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ric objects (stars) living in a space filled with the Chaplygin gas. Two cases are considered. In 
the normal case all solutions (excluding the de Sitter one) realize a three-dimensional spheroidal 
geometry because the radial coordinate achieves a maximal value (the "equator"). After crossing 
the equator, three scenarios are possible: a closed spheroid having a Schwarzschild-type singularity 
with infinite blue-shift at the "south pole", a regular spheroid, and a truncated spheroid having a 
scalar curvature singularity at a finite value of the radial coordinate. The second case arises when 
the modulus of the pressure exceeds the energy density (the phantom Chaplygin gas). There is no 
more equator and all solutions have the geometry of a truncated spheroid with the same type of 
singularity. We consider also static spherically symmetric configurations existing in a universe filled 
with the phantom Chaplygin gas only. In this case two classes of solutions exist: truncated spheroids 
and solutions of the wormhole type with a throat. However, the latter are not asymptotically flat 
and possess curvature singularities at finite values of the radial coordinate. Thus, they may not be 
used as models of observable compact astrophysical objects. 
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I. INTRODUCTION 

Because of the nowadays accepted existence of cosmic 
acceleration [l|, the study of spherically symmetric 
solutions of the Einstein equations [3, [J] in presence of 
dark energy is of much interest. This study has already 
been undertaken for instance in H, H, H, Q . 

One of the simp lest models for dark energy is the Chap- 
lygin gas [l(| Elf- The model is based on a perfect fluid 
satisfying the equation of state p = —A/p, where p is 
the pressure, p is the energy density and A is a positive 
constant. Some studies have already appeared where the 
problem of finding spherically symmetric or wormhole- 
like solutions of Einstein's equations in presence of the 
Chaplygin gas have been addressed d, 0] ■ 

Here we study static solutions of the Tolman- 
Oppenheimer-Volkoff (TOV) equations for spherically 
symmetric objects living in a space filled with the Chap- 
lygin gas. Results obtained appear to be very different 
from the apparently similar problem of stars in presence 
of a cosmological constant Q . Indeed, in the latter case 
the exterior solution of the TOV equations is nothing but 
the well-known Schwarzschild-de Sitter geometry, while 
the interior problem essentially coincides with the stan- 



dard TOV case. The only difference is that the pressure 
does not vanish at the star surface; on the contrary, it is 
negative and its absolute value is equal to the cosmolog- 
ical constant. 

Instead, the Chaplygin gas strongly feels the presence 
of the star, and consequently the solution acquires quite 
unusual features. These features are the existence of a 
maximal value of the radial coordinate, dubbed "equa- 
tor" , and the appearance of curvature singularities at 
some finite values of the radial coordinate r. We find 
also that for the case of the phantom Chaplygin gas when 
the absolute value of the pressure is greater than the 
energy density (|p| > p) wormhole-like solutions with a 
throat exist. However, these solutions cannot be identi- 
fied with the usual Morris-Thorne-Yurtsever wormholcs 
[l2| because they are not asymptotically flat. Moreover, 
they possess curvature singularities at finite values of r 
as well. 

The structure of the paper is as follows: in Sec. II 
we write down the TOV equations in presence of the 
Chaplygin gas and describe two exact solutions of them; 
Sec. Ill is devoted to the analysis of the normal case 
where \p\ < p while in Sec. IV we study the phantom case 
with \p\ > p, and we consider the star-like configurations 
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analogous to those studied in Sec. Ill; in Sec. V we 
consider the solutions of the TOV equations existing in 
a universe filled exclusively with the phantom Chaplygin 
gas. Sec. 6 contains conclusions and discussion. 



II. TOLMAN-OPPENHEIMER-VOLKOFF 
EQUATIONS IN THE PRESENCE OF THE 
CHAPLYGIN GAS 

We suppose that the universe is filled with a perfect 
fluid with energy-momentum T^ u = (p + p)u^u v — g^ v p 
and consider a static spherically symmetric interval 

ds 2 = e^dt 2 - e^dr 2 - r 2 (d9 2 + sin 2 6dc/) 2 ). (1) 

Then the Einstein system reduces to the following pair 
of equations: 



(2) 
(3) 

(4) 
= 1 
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8irp, 

plus the energy-momentum conservation equation 

dp dv p + p 
dr dr 2 

Solving Eq. ^ with the boundary condition e _Al (°) 
gives 



where, as usually, M(r) = Air drr 2 p(r 



alent to 



dM , , , 

— = 4nr 2 p, M(0) 
dr 



0. 



(5) 

This is equiv- 
(6) 



Eq. ©, @ and_|5|_together give rise to the TOV differ- 
ential equation 



dp (p + p) (M + 4irr 3 p) 
dr ~ r(r- 2M) ' 



(7) 



Complementing equations ([6]) and ([7]) with an equation 
of state relating p and p one has a closed system of three 
equations for the three variables p, p and M. In this 
paper we investigate the case when the fluid is the Chap- 
lygin gas whose equation of state is 



P 



A 2 
P ' 



(8) 



Then Eqs. © and ([7]) give rise to the following system 
of first-order differential equations for p and M: 



dp 
dr 
dM 
dr 



(A 2 — p 2 )(M + 47rr 3 p) 



pr(r 
4ttAV 
P 



2M) 



(9) 
(10) 



We denote the radius of the star by rt,. As usual we 
suppose that the pressure is continuous at the surface of 
the star. The "exterior" problem amounts to considering 
system (|10p in the interval r > r& with some properly 
chosen boundary conditions and M(r&) at r = r&. 

It is easy to see that at r > r^ the system admits two 
exact solutions with constant pressure. The first of them 



p = -p = -A, M = -7rAr 3 , 
3 

describes the geometry of the de Sitter space with 



e v = e - " = 1 



' dS 



8ttA 



The second solution is the Einstein static universe 
A 

with the radius 



V3A, M = 



'VI 

8ttA 



rds 
3V4- 



(11) 



(12) 



(13) 



(14) 



III. THE NORMAL CASE: |p] < p 

We now consider solutions with a non-constant pres- 
sure. Some additional constraints have to be imposed on 
the boundary conditions: 



A < p{r b ) < 0, 



M{ n ) < 



(15) 



(16) 



First of all, note that the pressure p cannot attain the 
values p = and p = —A in the region where 2M(r) < 
r. Indeed, for 2M(r) < r the right-hand side of ([9]) is 
negative, while in order to approach p — starting from 
negative values of p, it is necessary to have dp/ dr > 0. 
In addition, let us rewrite Eq. © as follows: 



dln(A 2 -p 2 



-2dr 



(M + Airr^p) 
r{r-2M) 



(17) 



and suppose that p — > — A as r — > n, with r% > 2M[r\). 
Then, upon integration of Eq. (TTT|) , we get a divergence 
at the left-hand side and a regular expression at the right- 
hand side, a contradiction. 

Thus, as long as the condition 2M(r) < r is satisfied 
we have —A < p(r) < and p(r) > A. This means that 
the mass M(r) is growing at least as fast as r 3 , so that 
at some radius r — ro the equality 



M(r ) 



111 
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(18) 
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is achieved. Then at r = ro we must have p(ro) = po = 
— 1 5 . Indeed, let us expand the relevant quantities 



around rp: 



M(r) 



- 

2 



M(e) 



p(r) = po + P», 

where M(e) and p(e) tend to zero when e - 
© has the following asymptotic form: 

dp = (A 2 -p 2 )(l + 8nr 2 p ) 
de 2e(p + 8^ A 2 ) ' 



(19) 
(20) 
(21) 

0. Equation 
(22) 



from which it is easily seen that p ~ lne when e — > 0, 
unless po = —3-^-7- 

We now determine a lower and an upper bound for the 
radius tq. First, note that since po > —A we have 



''0 > 



1 



8ttA 



(23) 



On the other hand, since p(r) > A, Eq. ([6]) implies that 
of M{r) > so that 



8ttA 



The asymptotic equation for p has the form 
Ts-Ye +Co 

where 

C a = 

and its solution is 



1 

8^ 



32tt 2 A 2 ^ 



p = Ay/e + 2C e 



(24) 



(25) 



(26) 



(27) 



where A is an arbitrary coefficient. Thus, the family of 
the solutions p(r), M(r) can be characterized by the two 
parameters ro and A which, in turn, are determined by 
the boundary conditions M(rb),p(rb) on the surface of 
the star. 

However, the coordinates which we have used so far 
are not convenient for the problem under consideration 
since the metric coefficient g rr = = (1 — 2M/r)~ 1 has 
a fictitious (coordinate) singularity at r — r§. Therefore, 
instead of the coordinate r, we introduce a new coordi- 
nate x defined by r = ro sinx, so that the corresponding 
metric becomes 



ds 2 



^dt 2 



e^ x >dx z - r 2 sin 2 X (d0 2 + sin 2 



(28) 



Then the tt-component of the Einstein equations has 
the form 

e' /i (/2 / cotx + 2-cot 2 x)+ 9 1 2 =8tt P (29) 

rgsm x 

where prime denotes differentiation with respect to the 
variable x- Integration with the boundary condition 
e -p(o) _ \/r 2 gives 



1 



1 



87rr 2 f x 



r 2 cos 2 x \ sin x „/ 



P(x) sin 2 xcosxrfx 



From Eqs. (|T8J and it follows that 

r 2 = 1 

8lr Jo Pbti sin2 * cos 



(30) 
(31) 



that implies the positivity and finiteness of the expression 

dsn}. 

For latter purposes we also write down the XX compo- 
nent of the Einstein equations: 



= A (cot 2 x + ^'cot%) 



1 



r 2 sin 2 x 



871-p. (32) 



The energy-momentum conservation equation now 
reads 



2p' 
p + p' 



(33) 



and for the case of the Chaplygin gas it can be easily 
integrated to give 



A 



A 2 -p 2 



(34) 



where A is some positive constant which fixes the choice 
of the time scale. 

Combining Eqs. ([29]), (J32]) and ((33]) one gets 



where 



p , = cos x(p + p) (M + 4?rrg sin 3 X p) ^ 



sinx(ro sinx — 2M) 



Af(x) = 47rrg/ p sin 2 X cos X dx (36) 
Jo 



and 



W = Airrfp sin 2 x cos x, M (0) = . (37) 

Relation |[T8"]) can be rewritten as 

r = 2M(tt/2), (38) 

and for the Chaplygin gas Eqs. ([33]) and ([57]) acquire the 
forms 

/ _ cos x(A 2 - p 2 )(M + 47rr 3 sin 3 X P) ,„ Q v 

P ~ '• 1 '• nl7\ ' V > 

p sm x( r o sm X ~ ) 
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M' 



47rrQA 2 sin xcosx 



P 



(40) 



respectively. 

We study these equations in the vicinity of the "equa- 
tor" (x = 7r/2) by introducing a small positive variable 
a such that 



7T 



and the functions p(a) and M(a): 

1 



(41) 



8^r 2 



M 



ro 
2 



M(a). 

A simple calculation shows that 

M = lQn 2 rlk 2 a 2 + • 



(43) 



(44) 



while for p(a) one can write down the following asymp- 
totic equation 



dp p 
da a 



Cia, 



where 



| - 327T 2 7-4A 2 

8^ ' 



The solution of Eq. (@5j) is 

p = Ba + C\a 2 , 



(45) 



(46) 



(47) 



with B an arbitrary constant. This solution can be con- 
tinued to negative values of the parameter a that cor- 
responds to the equator crossing. Thus, all trajectories 
intersecting the equator x = 7r/2 can be characterized 
by two parameters, which could be chosen as r and B. 
The static Einstein solution corresponds to the values 
B = and C\ = 0, the latter condition being equivalent 
to r 4 = 3/64tt 2 A 2 . 

In order to investigate the behavior of the trajecto- 
ries after crossing the equator, we find it convenient to 
introduce a new variable 



y 



l 



smx 



(48) 



so that 1 < y < oo. In terms of this variable Eqs. 
([3i||), (HOI can be rewritten as 



dp 
dy 



(A 2 — p 2 )(My 3 + 47rrgp) 



py 3 (r - 2My) 



dM 
dy 



47rA 2 r-3 
py 4 



(49) 



(50) 



Now one can show that the expression ro — 2My in the 
denominator of the right-hand side of Eq. (|4"9"1) is always 
positive at y > 1. In order to prove this statement, we 
first show that it is true if p > —A. 

To this end we introduce the function 



f(y) =r - r y + 



8^Arg(y 3 -l) 
3y 2 



Since p > —A, it satisfies the inequality 



(42) We have 



f(y) <r - 2M(y)y, 



/(1) = 



and 



so that 



f'(y) 



8?rArg 167rArg 



> 0. 



(51) 



(52) 



(53) 



(54) 



(55) 



Now assume that the function r —2M(y)y becomes equal 
to zero at some y = y\ > 1. This means that at some 
value y = yi < yi, the function f(y) vanishes. In turn, 
this last condition requires the vanishing of the derivative 
f'(y) at some value y = y% <yi- From Eq. (|5l)l one finds 



2/3 



167rAr 



3 - 87rAr 2 



1/3 



and the requirement y^ > 1 is equivalent to 



1 < 87TA7- 2 , < 3. 



(56) 



(57) 



The vanishing of the function tq — 2My at the point y\ 
implies also the vanishing of the expression My 3 + Anr^p 
at this point, i.e. the vanishing of the numerator of the 
expression in the right-hand side of Eq. ([49|) . Thus, the 
values of the functions M and p at the point y\ are given 
by 



M(yi) 



TO 

2yi 



p(yi) 



StTT 2 ' 



It follows from the condition p > — A that 
y\ < 87rAr 2 , 

while from 1/3 < y± we find 

y\ < SttAt- 2 . 



(58) 



(59) 



(60) 



Substituting the expression (|56|) into the inequality ([601 
we get that this inequality is satisfied provided 



87rAr 2 > 4 



(61) 
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that contradicts the condition (|57|) . 

Finally, looking at Eq. one can see that the 

pressure can, in principle, achieve the value p = —A 
at some value y = y\ only if ro/y\ = 2M(y\) and 

{j—jr 2 - — lj > 0. A simple analysis similar to the one 

carried out above shows that this is impossible as well. 
Thus, we have shown that the expression ro — 2M(y)y 
cannot vanish at any value of y in the range 1 < y < oo. 

We now study the behavior of the pressure at y > 1. 
Here we find three families of solutions (geometries) . The 
first one contains trajectories arriving at the south pole 
of the three-dimensional spatial manifold (y = oo, \ — 51 ) 
with some value > p(oo) > —A. Looking at Eq. (|4"9")) 
we see that the necessary condition for such solutions to 
exist is the convergent behavior of the integral 



dy 



My 3 + Airrlp 
y 3 (r ~2My) 



(62) 



at y — > oo that implies the vanishing of the function M 
at y — ► oo. Indeed, assume M(oo) = Mq ^ 0. Then 
Mo > contradicts the positivity of the expression r — 
2My, while M < implies the integral (f6"2"|) to diverge 
logarithmically. Thus, the only value M(oo) compatible 
with p(oo) > —A is M (oo) = 0. Then it follows from Eq. 
(|50f that the asymptotic behavior of M at y — > oo is 



m 

M = -s, 

yA 



where 



47rA 



3p(oo) 



(63) 



(64) 



Substituting the value of m into the integrand of the 
right-hand side of Eq. , we see that the sign of the 
derivative dp/dy at y — > oo is determined by the sign 
of the expression (A 2 — 3p 2 (oo)). If p < — A/V3 this 
derivative is negative, while it is positive for p > —A/ 

Thus, there exists a two-parameter family of regular 
spacetime geometries for which the spatial manifold rep- 
resents a three-dimensional spheroid parameterized by 
the two parameters ro and p(oo) and the metric coeffi- 
cient g u given by the formula (|34|) is always positive. 

The second family of geometries includes those where 
the value of the pressure p becomes equal to zero. Let us 
describe basic features of such geometries. We suppose 
that p(yo) = at yo > 1- The function M(y) cannot be- 
come negative at y = yo because in this case the deriva- 
tive dp/dy would be negative and it would be impossible 
to reach the value p(yo) = 0. Hence we consider the case 
when M(yo) = M where < M < £k We assume 
that in the neighborhood of the point yo the pressure 
behaves as 



where D and a are some positive constants. Substituting 
the expression (|65|) into Eq. (|49|) . one gets 



D 



I 2A 2 M 
r Q - 2M y ' 



(66) 



Consider also the case when M(yo) = 0. In this case we 
look for the expressions describing the behavior of p and 
M in the vicinity of y — yo in the form 



M(y)=M 1 (yo- y y 



p(y) = -E(vo - y) 1 -, 



(67) 



(68) 



where Mi and E are positive and < (3 < 7. Substi- 
tuting expressions (|6"T)) and (|6"8|) into Eqs. (gH), ([50)1 one 
finds the following values for the parameters /?, 7, Mi and 
E: 



E : 



^ 187rA 4 r 2 



yt 



1/3 



Mi 



2r 
3A 2 



E- 



(69) 

Note that in this case the values of the of the parame- 
ters Mi and E are uniquely fixed by the value of yo- In 
the case of M(y ) > considered above one has a one pa- 
rameter family of geometries parameterized by the value 
of Mq or by the value of D. Thus, it seems that one 
has a three-parameter family of geometries correspond- 
ing to p — * and these parameters are ro,yo and Mq, 
However, Eqs. ((66|) or ([69]) describe necessary conditions 
which should be satisfied to provide the existence of the 
geometry having the maximal radius ro and the pressure 
p vanishing at y = yo- Not all the solutions satisfying the 
relations (|66|) or (|69|) correspond to "initial conditions" 
at y = 1, i.e. M(y = 1) = ro/2. Moreover, taking into 
account the monotonic behavior of the function M(y) 
one can believe that at least one value of the parameter 
M(yo) corresponds to a geometry with the desirable ini- 
tial and final conditions. Thus, the family of solutions 
(geometries) ending with p — is also two-parametric 
and can be parameterized by the two parameters ro and 

yo- 

An interesting feature of the geometries described 
above consists in the presence of the singularity at y = yo . 
Indeed, the Chaplygin gas equation of state implies an 
infinite growth of the energy density when the pressure 
tends to zero that in turn determines the divergence of 
the scalar curvature. Thus, the spacetime under consid- 
eration cannot be continued beyond y — yo or, in other 
terms, beyond x = t — arcsinj/Q -1 . 

The third family of possible geometries includes those 
for which the pressure p tends to the value —A when 
y — > 00 (x — > 7r). In this case the acceptable behavior of 
the function M is M(oo) = -M 2 , where M 2 > 0. The 
behavior of the pressure at y — > 00 can be represented as 



p(y) = -D(y - y) a (65) 



p(y) = -A + p(y) 



(70) 
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where p{y) is a positive function vanishing at y — > oo. 
Substituting ([70)) into Eq. (J35J), one gets an asymptotic 
equation 



dp 
dy 



P 



which solution is 



F 

P = — 

y 



(71) 



(72) 



where F is a positive constant. This constant does not 
depend on the value of the parameter M%. Thus, the fam- 
ily of the geometries with p — > —A at y — > oo appears to 
be described by the three parameters fo, -F and M2. How- 
ever, as in the case of the family of geometries described 
above with the pressure vanishing at some of y = yo, we 
are not free in the choice of the value of M2 after the 
values ro and F are fixed. Indeed, due to the monotonic 
behavior of the function M(y), at least one value of the 
parameter M2 corresponds to a solution M(y) satisfy- 
ing the initial condition M(l) = r /2. Thus, one has a 
two-parameter family of geometries defined by fixing thee 
values of ro and F. These geometries have a singularity 
of the Schwarzschild type at y = 00 (x = 7T, r = 0) due 
to the nonvanishing mass M = —Mi- They have another 
curious feature: the metric coefficient gu = e vtyX ^ given 
by the formula (|34[) tends to infinity as p — > —A and, 
hence, intervals of the proper time dr — y/~g~[idt tend to 
infinity. So, one has an infinite blue shift effect in con- 
trast to the well-known red shift effects in the vicinity of 
the Schwarzschild and de Sitter horizons. 

Summarizing, solutions of the Tolman-Oppenheimer- 
Volkoff equations in the presence of the Chaplygin gas 
have the following curious features: 

1. All the spatial sections of the spacetime manifolds 
(excluding a special case of the de Sitter spacetime) 
are closed. 

2. Some geometries have a divergent scalar curvature 
invariant at a finite value of r. 

3. Some geometries manifest an infinite blue-shift ef- 
fect. 

Unfortunately, the relations between the boundary 
conditions p b ,M b , the parameters, characterizing the 
crossing of the equator ro, B and the "final parameters" 
characterizing the three family of geometries with quali- 
tatively different behaviors at \ > 7r /2 cannot be found 
analytically and should be studied numerically. 



IV. THE PHANTOM CASE: \p\ > p 

Now consider the system of equations ©, (|10j) with 
the boundary condition 



p{n) < -A. 



(73) 



In this case \p\ > p that corresponds to phantom dark 
energy and, in principle, to a possibility of creation of 
wormholcs. If the condition ([75)) is satisfied, two cases 
are possible. 
Case A: 



M(r b ) + Anrlp{r b ) < 0. 



(74) 



In this case, the pressure p is decreasing and its absolute 
value is growing. Correspondingly, the energy density 
is also decreasing and hence the mass M grows slower 
than r 3 . Then the left-hand side of the expression ([74]) is 
decreasing, too, so the expressions (M + 4irr 3 p), (r — 2M) 
and (p 2 — A 2 ) cannot change their signs. 
We examine three possible subcases: 

1) p tends to some finite value —00 < p\ < —A when 
r — > 00; 

2) p tends to —00 when r — > 00; 

3) p grows indefinitely when r tends to some finite value 
r\- 

The subcase 1 cannot take place because the left-hand 
side of Eq. © is regular while its right-hand side diverges 
as r 2 when r — > 00. 

Likewise, the subcase 2 cannot be realized as well. 
Indeed, suppose that p — —p\r a ,a > 0,pi > when 
r — > 00. Then Eq. ([74")) becomes 



dp 
dr 



-47rp 2 r, 



which implies 



a - 1 = 2a + 1 



(75) 



(76) 



or a = — 2 which contradicts to the positivity of a. 

We are left with the subcase 3 which can be realized 
with 

Pi 



P : 



7*1 — r 



-, pi > 0, when r — ► r\. 



(77) 



Substituting expression ([77)) into Eq. ([9]), we have the 
following relation: 



Pi 



pi 



{n - R) 2 {n - r) 2 n - 2M l 



(78) 



from which we find the value of parameter pi as a func- 
tion of the radius 7*1 and the mass M(ri) = Mi: 



Pi 



n - 2M1 

47rr 2 



(79) 



Thus, we obtain a two-parameter family of solutions in 
which one encounters a singularity at r = r\ because the 
scalar curvature R diverges there. 
Case B: 



M(r b ) + ^r 3 b p(r b ) > 0. 



(80) 



In this case the pressure grows with r. Then, in the right- 
hand side of Eq. © we have three decreasing positive 
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terms (p 2 -A 2 ), (M + Anpr 3 ) and (r — 2M). The problem 
is which of them vanishes before the others, if any. 

The above terms cannot simultaneously remain posi- 
tive as r — > oo because in this case \p\ > A > p and the 
expression (M + Anpr 3 ) would unavoidably change sign. 

The case when p = —A while the other two expressions 
remain positive is also excluded. Indeed, if p — ► —A as 
r approaches some finite value, the left-hand side of Eq. 
(jTTJ) has a logarithmic divergence while its right-hand 
side is regular. On the other hand, if p — > —A as r — > oo, 
the expression (M + iwpr 3 ) will change its sign. 

If (p 2 — A 2 ) and (r — 2M) vanish at some r = ro, Eq. 
(fT?"]) takes the form 



dln(p 2 - A 2 ) = -- 



dr 



ro 



which implies 



P 



A 2 



r -r 



(81) 



(82) 



contradicting the hypothesis. Thus, the pressure cannot 
achieve the value p = — A. 

As shown in the preceding section, the denominator 
(r — 2M) at the right-hand side of Eq. © can only 
vanish at some r = ro simultaneously with (M + Airpr 3 ) 
and the pressure at r = r should be equal p = —%^pz 

(equator). Let us prove that it is impossible to achieve 
the equator if p < —A. Indeed, the formula (jTQJ) shows 
that in this case 

r - 2M = r + (r - r ) - 2M(r ) - 2M'(r ) (r - r ) 
= (r-r )(l- 647r 2 A 2 r4) (83) 

as r — » ro. The difference (r — 2M) should be positive 
as r — * r from below, so the expression (1 — 647r 2 A 2 rQ) 
should be negative. However, at \p\ > A this expression 
is positive - a contradiction. 

Thus, if (|80p is satisfied, p grows until some maxi- 
mum value p m ax < —A when the expression (M + 4irpr 3 ) 
changes sign while the terms (p 2 — A 2 ) and (r — 2M) are 
positive, and we led back to the case A. 

So, we have proved that only two regimes are pos- 
sible for a star-like object immersed into the phantom 
Chaplygin gas. If initial conditions satisfy ([74")) , the pres- 
sure is decreasing and diverges at some finite value of 
r. Then the space-time acquires a scalar curvature sin- 
gularity there. On the other hand, if initial conditions 
satisfy ([50)) . then the pressure grows with the r until 
some maximum value p = p ma x where the expression in 
the right-hand side of Eq. ([9]) changes sign. After that 
we come back to case A: the pressure decreases and ex- 
plodes according to ([77)) at some finite value of the radial 
coordinate r. No equator and no horizon are attained 
in the case of the phantom Chaplygin gas: the quantity 
(r — 2M) is always positive. As in the case of the non- 
phantom Chaplygin gas, the relation between the initial 
values of the parameters functions p(rb) and M(rb) and 



the parameters r\ and M(n) cannot be found analyti- 
cally. 

We can summarize the results of the above considera- 
tions in the following 

THEOREM. In a static spherically symmetric uni- 
verse filled with the phantom Chaplygin gas, the scalar 
curvature becomes singular at some finite value of the 
radial coordinate and the universe is not asymptotically 
flat. 



V. STATIC SPHERICALLY SYMMETRIC 
UNIVERSE FILLED EXCLUSIVELY WITH THE 
PHANTOM CHAPLYGIN GAS 

Now we study spherically symmetric static solutions 
for a universe filled exclusively with the phantom Chap- 
lygin gas. The theorem above is valid in this case, too. 
This situation is of much interest because when the weak 
energy condition is violated, p + p < 0, wormholes may 
appear (though not necessarily, see 0] in this connec- 
tion). Suppose that at some finite value of the radial 
variable r = r/,, the factor r^ — 2M(rb) is positive and 
p( r b) < — A. Then, as in the preceding analysis, one can 
consider the evolution of the functions M(r) and p(r) 
in accordance with Eqs. ([9]), (|10p but with a decreasing 
value of the radial variable r. Now only two possibilities 
may be realized: one can arrive at the value r = keep- 
ing always a positive value of the factor r — 2M, or one 
can encounter a situation when at some finite value of 
r = ro this factor vanishes. 

Consider first the case when (r — 2M) is positive for 
all values r > 0. Here, one can imagine two differ- 
ent regimes as r approaches zero. In the first one the 
mass in the vicinity of r = is positive and behaves as 
M ~ r a ,a > 0. In the second regime the mass tends 
to a negative constant when r — > 0. A detailed analysis 
shows that only the first regime is compatible with the 
TOV equations ([5]), (jTTJj) in the presence of the phantom 
Chaplygin gas. Precisely, in the vicinity of r = the 
pressure and the mass functions have the form: 



P = Po 



8^ 2 (3p 2 -A 2 )(p 2 -A 2 )r 2 



M 



4-7rA 2 r 3 
3po 



(84) 



(85) 



where po is an arbitrary number such that po < —A. Be- 
ing regular at the center r = 0, this static configuration 
develops a singularity at some finite value of the radius 
ri, where the pressure becomes equal to minus infinity. 
Thus, we have a one-parameter family of static spheri- 
cally symmetric solutions of the Tolman-Oppenheimer- 
Volkoff equations in a world filled with the phantom 
Chaplygin gas. This family is parameterized by the value 
of the pressure at the center r = 0. 
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Now suppose that the factor r > 2M(r) vanishes at 
some value r = r . An analysis similar to the one car- 
ried out in Sec. 3 shows that this is positive only if also 
the expression M + 4irpr 3 in the numerator of the right- 
hand-side of Eq. ([9]) vanishes at ro. In Sec. 3 the surface 
r = ro was called equator because it corresponded to the 
maximal value of the radial variable r. Now it corre- 
sponds to the minimal value of r, and it is nothing but a 
throat. Just like in the case of the equator considered in 
Sec. 3, the throat can be achieved only at p = — o-^t. In 
the phantom case p < — A, hence, there is a restriction 
on the size of the throat 



ro < 



1 



8ttA 



(86) 



In order to describe the crossing of the throat, it is con- 
venient to introduce the hyperbolic coordinate r\ instead 
of the radius r, 



ro cosh 7/, 



(87) 



which plays a role similar to that played by the trigono- 
metrical angle x m the description of the equator. 
Now the TOV equations look like 



dp (A 2 - p 2 )(M + iirpr^ cosh 3 77) sinh 
drj p cosh 77(7-0 cosh 77 — 2M) 

dM 47rA 2 rQ cosh 2 77 sinh 77 
drj p 



V 



(89) 



The solution of (|89[) at small values of 77 is 

M = r + 167r 2 A 2 ^?7 2 . (90) 
Representing the pressure as 

1 



p(rj) 



87T7- 2 



(91) 



one can rewrite Eq. (|55|) in the neighborhood of the 
throat as 



— = - + C T rj 
drj rj 

where the negative constant Ct is equal to 
647r 2 A 2 r 1 - 3 



Cn 



167T7- 2 



The solution of Eq. ^ is 



p = Drj + -Crrf 



(92) 



(93) 



(94) 



where D is an arbitrary constant. 

Thus, geometries with a throat constitute a two- 
parameter family, characterized by the value of the throat 
radius ro and by the value of the parameter D. From Eq. 



(194)1 we see that when the hyperbolic parameter 77 grows, 
the negative term ^CtV 2 starts dominating, while the 
pressure decreases and achieves an infinite negative value 
at a finite value of the radius r where we encounter a cur- 
vature singularity. The peculiarity of this configuration 
with a throat consists in the fact that these singulari- 
ties are achieved at both sides of the throat or, in other 
words, at one positive 771 > and one negative 772 < 
values of the hyperbolic parameter. The values 771 and 772 
correspond to values r\ = r cosh 77! and r 2 = r cosh 772 
of the radial parameter. If D = one has 7*1 = r 2 . The 
solutions with a throat could be seen as wormholc-likc, 
but in contrast with the Morris-Thorne-Yurtsever worm- 
holes [12], they do not connect two asymptotically flat 
regions. 

We summarize the preceding considerations concern- 
ing solutions of the Tolman-Oppenhcimer-Volkoff equa- 
tions in the universe filled with the phantom Chaplygin 
gas as follows. Suppose we start from an initial condition 
f—2M{¥) > and p(f) < —A for some given value f of r. 
Then letting the functions M(r) and p(r) evolve to val- 
ues r > f, we unavoidably arrive to an infinite negative 
value of the pressure at some finite value r / of the radius 
r, thus encountering a singularity. Instead, for r < f two 
qualitatively different situations may arise: either we ar- 
rive to r = in a regular way, or we may discover a throat 
at some finite value r = ro, this being the generic situ- 
ation. Upon passing the throat, an observer finds itself 
in another patch of the world and then, with the radius 
increasing, stumbles again upon a curvature singularity 
at some finite distance from the throat. 

Note that in the traditional view of wormholes, one 
supposes that there is a minimal value of the radial pa- 
rameter characterizing a throat, and the space-time at 
both sides of the throat is either asymptotically flat, or 
has some other traditional structure (for example, worm- 
holes could also connect two expanding asymptotically 
Friedmann universes). Here we have found a different 
kind of wormhole-like solutions: those connecting two 
space-time patches which have a scalar curvature singu- 
larity at some finite value of r. 

We conclude this section with the brief comment on the 
results of paper [7;] where the problem of the existence of 
wormhole solutions supported by the phantom Chaply- 
gin gas was studied. The main part of this paper is de- 
voted to the consideration of the so called "anisotropic" 
Chaplygin gas, i.e. a fluid whose radial pressure satisfies 
the Chaplygin gas equation of state, while the tangen- 
tial pressure can be arbitrary. Then the system of TOV 
equations is under-determined and its solution contains 
one arbitrary function. Choosing this function in a con- 
venient way, one can construct a lot of solutions, satisfy- 
ing the desired properties. However, this fluid is not the 
Chaplygin gas and, moreover, is not a barotropic fluid. 
The case of the isotropic Chaplygin gas is also considered 
in Ref.Q. The author studies numerically the behavior 
of metric coefficients in the vicinity of the throat without 
considering the problem of continuation of this solution 
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to larger values of the radial coordinate. However, the 
general theorem proved at the end of the section IV of 
the present paper states that all spherically-symmetric 
solutions supported by the phantom Chaplygin gas (with 
or without a throat) have a curvature singularity at some 
finite value of the radial coordinate and hence cannot be 
asymptotically flat. 

VI. CONCLUSION 

In this paper we have studied the Tolman- 
Oppenheimer-Volkoff equations for static spherically 
symmetric objects immersed in the space filled with the 
Chaplygin gas. Both cases, phantom and non-phantom, 
were considered. In the non-phantom case all solutions 
(excluding the de Sitter one) represent a spheroidal ge- 
ometry, where the radial coordinate achieves a maximal 
value (equator). After crossing the equator, depending 
on the boundary conditions, three types of solutions can 
arise: a closed spheroid having a Schwarzschild-type sin- 
gularity with infinite blue-shift at the "south pole" , a reg- 
ular spheroid, and a truncated spheroid having a scalar 
curvature singularity at a finite value of the radial coor- 
dinate. 

For the case of the phantom Chaplygin gas, the equator 
is absent and all star-like external solutions have the ge- 
ometry of a truncated spheroid having a scalar curvature 
singularity at some finite value of the radial coordinate. 
Besides, we have also considered the static spherically 
symmetric configurations existing in a universe filled ex- 



clusively with the phantom Chaplygin gas. Here two 
cases are possible: geometries which are regular at the 
center r = and having a scalar curvature singularity at 
some finite value of r, and geometries containing throats 
connecting two patches of the world which again have 
scalar singularities at some finite values of the radius. 
Because of these singularities, the construction of stable, 
traversable and asymptotically flat wormholes using the 
phantom Chaplygin gas is prohibited, in spite of breaking 
the weak energy condition in this case. 

Finally, we note that many of the solutions of the TOV 
equations studied above possess singularities arising at 
some finite values of the radial coordinate, while intu- 
itively it is more habitual to think about singularities 
arising at the point characterized by the vanishing of 
this coordinate. Something similar happens in the study 
of isotropic and homogeneous cosmological models, too. 
Here, in addition to the traditional Big Bang and Big 
Crunch singularities, an intensive study of the singulari- 
ties which take place at finite or at infinite values of the 
cosmological scale factor is under way (see e.g. [Hj|). 
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